Fermi surface topology of LaFePO and LiFeP 
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We perform charge self-consistent LDA+DMFT (density functional theory combined with dy- 
namical mean field theory) calculations to study correlation effects on the Fermi surfaces of the 
iron pnictide superconductors LaFePO and LiFeP. We find a distinctive change in the topology of 
the Fermi surface in both compounds where a hole pocket with Fe d z 2 orbital character changes its 
geometry from a closed shape in LDA to an open shape upon inclusion of correlations. The opening 
of the pocket occurs in the vicinity of the F" (Z) point in LaFePO (LiFeP). We discuss the relevance 
of these findings for the low superconducting transition temperature and the nodal gap observed in 
these materials. 
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The iron pnictides are an intensively studied new class 
of superconductors with superconducting transition tem- 
peratures up to currently 56 KP They show a variety 
of distinct properties in the normal as well as supercon- 
ducting state like absence/presence of magnetic ordering, 
weak to strong electronic correlations and nodal/nodeless 
superconducting gaps. Most parent compounds exhibit 
a non-superconducting ground state with antiferromag- 
netic order and become superconducting upon doping or 
application of external pressure. 

LaFePO is the iron pnictide in which superconductivity 
was reportecP for the first time at a critical temperature 
of about 6 K. In LaFePO superconductivity arises with- 
out doping or application of pressure, there is no long- 
range magnetic orderj^] and the superconducting gap is 
nodalfSHZl As a measure of electronic correlations, mass 
enhancement values have been reported from a number 
of different experiments including angle-resolved photoe- 
mission spectroscopy (ARPES)P optical conductivity, 9 
de Haas- van Alphen (dHvA)E3 and specific heat mea- 
surements}^ all point to a range m* j ttilda ~ 1. 7 — 2.2 , 
consistent also with existing LDA+DMFT studies PEH 

Recently, LiFeAs received particular attention due to a 
number of features that make it very attractive for both 
theoretical and experimental studies: like LaFePO, it ex- 
hibits superconductivity without doping or pressure at 
T c = 18 K 14 and is nonmagnetic; furthermore, it cleaves 
between adjacent Li layers mitigating the issue of polar 
surfaces for surface- sensitive probes. 

LiFeP shares these properties, yet its lower supercon- 
ducting transition temperature of 5 rd^ didn't trigger 
as much research. LiFeP shows contrasting behavior to 
LiFeAs in some respects, though. Most importantly, its 
superconducting order parameter is nodaP^ compared to 
the nodeless gap in LiFeAs. This is unexpected as pre- 
vious works suggested a relation between the lack of a 
third hole pocket at the T point and the formation of gap 
nodes ESEHl b u t both LiFeAs and LiFeP show three hole 
sheets at V. Effective masses in LiFeP were extracted 
from resistivity and upper critical field measurement^^ 
where the mass enhancements are estimated to be smaller 



by a factor of ^2 compared to LiFeAs; this corresponds 
to a mass enhancement of 1.5-2 over the LDA value. 
De Haas-van Alphen experiments^! give values 1.6-3.3 
for the mass enhancements. However, so far no theo- 
retical studies investigating the effects of correlations on 
LiFeP have been reported. 

Thus, both compounds are considered rather weakly 
correlated. However, in the present work we argue that 
the inclusion of correlations has a profound impact on the 
Fermi surface topology of both materials which in the 
case of LaFePO agrees with ARPESS observations; for 
LiFeP, ARPES measurements are not yet reported and 
the available dHvA^data do not allow to unambiguously 
decide on the k z extension of the Fermi surface sheets as 
will be discussed below. The features presented here have 
not been touched upo n in t he reported LDA+DMFT cal- 
culations on LaFePOpMl for LiFeP we present, to our 
knowledge, the first LDA+DMFT study in the literature. 

We performed full charge self-consistent calculations 
following the scheme from Ref. [fusing a combination of 
electronic structure calculations in the full potential lin- 
earized augmented plane wave (FLAPW) framework as 
implemented in WIEN2k 24 with DMFT. For solving the 
impurity problem wc employed the hybridization expan- 
sion continuous-time quantum Monte Carlo meth ocP^ 
(CT-HYB) as implemented in the ALPS code}2£122l only 
density-density terms of the Hund's coupling were con- 
sidered. We performed calculations on all available ex- 
perimental crystal structures as reported in Refs. [3] and 
4 (LaFePO) and Refs. [151 and [22] (LiFeP) with space 
group P 4/nmm; results are shown for the structures 
from Ref. [3] (LaFePO) and Ref. HU (LiFeP). The energy 
window for the construction of a localized Wannier basis 
was chosen to range from -5.4 eV to 2.7 eV (-6 eV to 
3.15 eV) with respect to the Fermi energy for LaFePO 
(LiFeP). The Monte Carlo sampling was done at an in- 
verse temperature /3 = 40 eV -1 with 3 x 10 6 sweeps. 

For the interaction parameters, we use the definitions 
of U = F° and J = (F 2 + F 4 )/U in terms of Slater 
integrals^ F k with [7 = 4 eV, J = 0.8 eV, and the 
around-mean-fielcpSJ (AMF) double counting correction. 
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Some low energy features, in particular the size and shape 
of the Fermi hole pockets which are of central interest 
here, turn out to be rather sensitive to details of the 
LDA+DMFT calculation like the choice of double count- 
ing; however, we tested variations in these technical as- 
pects and found the reported features -while being af- 
fected quantitatively- qualitatively consistent among cal- 
culations with different double countings, interaction pa- 
rameters, and reported crystal structures. 

In the following, orbital characters are labeled in a co- 
ordinate system which is 45° rotated with respect to the 
crystallographic axes, i.e. x and y point to nearest Fe 
neighbors in the Fe-P plane. 

In Table U we list the orbital-resolved mass en- 
hancements for both compounds. The mass enhance- 
ments for LaFePO are in line with the measured val- 
ues ~1.7-2.2 from the various experiments^^ as well as 
previous LDA+DMFT st udies where mass enhancements 
~1.6-2.2 were calculatedPEH Note that for LaFePO the 
effective masses are higher for the e g orbitals whereas in 
LiFeP (and most other iron pnictides) the t<z g orbitals 
show stronger renormalization. This is because of the 
crystal field splitting which in LaFePO puts the ti g or- 
bitals below the e g orbitals thereby promoting a ground 
state with configuration e^t^g in the atomic limit. This 
suppresses inter-orbital fluctuations among the e stat es, 
rendering these orbitals more correlated in effectP^^ As 
a consequence, the d z i orbital is the most strongly cor- 
related one in LaFePO. 

The values for LiFeP range between 1.4 and 1.7 which 
is roughly a factor of 2 smaller than in LiFeAs where 
ARPES^ and dHvA^ experiments yield mass enhance- 
ments of 3-4. This is in agreement with resistivity and 
upper critical field measurements^ which also give a fac- 
tor 2 reduction with respect to LiFeAs. 



TABLE I. Orbital-resolved mass enhancements m'/wiLDA- 



Orbital 




d z 2 


^x 2 — y 2 


dxy 


^xz 1 yz 


LaFePO: 
LiFeP: 


m* 
m LDA 

m 
m LDA 


1.85 
1.52 


1.70 
1.39 


1.54 
1.71 


1.69 
1.62 



The momentum-resolved spectral function of the two 
materials is presented in Fig. [T] in comparison with the 
LDA band energies. The excitations around the Fermi 
energy are well-defined revealing the Fermi liquid nature 
in accordance with resistivity measurements!^ Most im- 
portantly, both compounds feature a band placed just 
below the Fermi energy in LDA which gets shifted above 
Ep upon inclusion of correlations in the vicinity of the 
r (Z) point in LaFePO (LiFeP). In both compounds this 
band has d z i orbital character around Ep for the path 
shown in Fig. [T] and originates from the hybridization 
with phosphorus p states. 

As shown in Figs. [2] and [3j this crossing of the Fermi 
energy is naturally accompanied by the appearance at 
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FIG. 1. (Color online) Momentum-resolved spectral function 
A(k, oS) together with LDA bands in the vicinity of the Fermi 
surface topology change. 



the Fermi surface of an outer hole pocket centered at L in 
LaFePO (see Fig.^b)) and an inner hole pocket centered 
at Z in LiFeP (see Fig. |3|b)). As has been suggested^ 
for the iron pnictides, appearance of a pocket with d z 2 
character may alter the superconducting pairing function 
to a nodal state and reduce the strength of the pairing as 
it is observed in LaFePO and LiFeP in contrast to their 
arsenic counterparts LaFeAsO and LiFeAs. 

Due to the hybridization with the phosphorus states, 
the position of the crossing band in LDA is sensitive to 
the phosphorus z position and doping. Therefore, al- 
though the large outer hole pocket appearing in LaF ePO 
has been observed in ARPE^P, it has been suspectecP^ 
that the appearance of this pocket is caused by surface 
doping indicated by a too small electron count obtained 
in ARPES. In contrast, our calculations yield the change 
in the Fermi surface topology as a result of electronic cor- 
relations only. Note that the total charge in the crystal 
is conserved in our calculations and the opening of the 
pocket (i.e. increase of the Fermi surface volume) merely 
amounts to a charge transfer from the d z i orbitals to the 
tig orbitals. Despite the sensitivity of the band posi- 
tion in LDA (for LaFePO the band energies differ by ap- 
prox. 12 meV between the two published structures^) 
we found the opening of the pocket in both structures and 
with very similar pocket sizes. As for the electron defi- 
ciency in ARPES, the huge size of the measured pocket 
(>12 kT as compared to <5 kT in our calculations, cf. 
orbit 3a in Fig. Qa)) probably still results from a charge 
effect. 

As a result, the calculated dHvA frequencies for 
LaFePO experience significant shifts as shown in 
Fig. RJa) . The frequencies correspond to extremal pocket 
sizes (orbits) that are observed at a given angle O with 
respect to the k z axis. The outer hole pocket experiences 
a large increase compared to LDA and the opening at 
Z adds a new frequency 3a for the minimal orbit. As a 
result of charge conservation, the outer electron pocket 
4a/b is also blown up. The enlargement of the electron 
pocket seen in our calculations is not observed in the 
dHvA experiment.^ The hole pocket itself is not mea- 
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FIG. 2. (Color online) Fermi surfaces of LaFePO in the k z = and k z — ir plane (left panels) and the k x — k y plane (right 
panels) for (a) LDA and (b) LDA+DMFT. The colors give the dominating orbital character. 
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FIG. 3. (Color online) Fermi surfaces of LiFeP in the k z = and k z = ir plane (left panels) and the k x — k y plane (rights 
panels) for (a) LDA and (b) LDA+DMFT. The colors give the dominating orbital character. The arrows indicate the fc-points 
for the calculation of the effective masses in Tab. [TX1 



sured in dHvA (7 out of the 10 predicted frequencies 
are present in the measurements). In LDA, inclusion 
of spin-orbit (SO) coupling reduces the size of the inner 
hole pocket la/b but SO coupling is not included in our 
LDA+DMFT calculations. It is therefore likely that this 
pocket shrinks even more than predicted by us, thereby 
reducing the total Fermi surface volume enclosed by the 
hole pockets; this could approximately compensate for 



the added volume from the opened hole pocket without 
enlargement of the electron pocket. 

For LiFeP, the inclusion of interactions induces only 
moderate changes in the sizes of the Fermi surface sheets, 
see Fig. Qb). The frequency shifts with respect to LDA 
are in qualitative agreement with the experimental dHvA 
data from Ref. the middle hole pocket 2a/b shifts 
most and shrinks in size by approx. 0.39 kT (compared 
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FIG. 4. (Color online) dHvA frequencies with respect to mag- 
netic field angle obtained within LDA (dashed thin lines) and 
LDA+DMFT (solid lines). Orbits 1, 2, and 3 refer to the 
inner, middle, and outer hole pocket, and orbits 4 and 5 to 
the outer and inner electron pocket. 



to approx. 0.95 kT in the experiment), the other sheets 
are less affected (orbit 4b also shrinks substantially but 
it it not measured in dHvA). The opening of the inner 
hole pocket in LDA+DMFT has almost no effect on the 
dHvA frequencies: also in LDA, two frequencies for the 
inner hole pocket are expected due to a weak peanut-like 
distortion which gives a minimal orbit around T and a 
maximal orbit at k z x /2. In LDA+DMFT, the size 
of the pocket at Z (the new maximal orbit) essentially 
equals the maximal orbit size from LDA, thereby mimick- 
ing the LDA orbits. Since the position of the orbits in the 
BZ cannot be determined from dHvA, no clear distinction 
between the predictions from LDA and LDA+DMFT can 
be made from the published data. Measurements up to 
O = 90° which in principle allow to differentiate between 
open and closed pockets (F cos 8 drops to zero for a closed 
pocket) have been performed and indicate a tendency of 
FcosO towards small values but are not conclusive be- 
cause of the very weak signalP^ 

While a good qualitative agreement is reached between 
our results and dHvA observations, our calculations do 
not lead to a sufficient shift to attain complete agreement 
with the experimental frequencies in LiFeP, in particu- 
lar the calculated reduction of the middle hole pocket is 
not pronounced enough. Note that spin-orbit coupling 
helps with the size reduction of this pocket, its effect is 
comparatively small here, though (about 0.2 kT). Limi- 
tations of our approach are also revealed by a comparison 
of the effective masses in LiFeP. The effective masses ob- 
tained from the dHvA measurements are rather uniform 
among all orbits except for the orbits 2a/b which show 
only half the mass enhancement of the other orbits. Since 



these mass enhancements refer to the Fermi surface or- 
bits rather than the localized orbitals, we calculated the 
LDA+DMFT mass enhancements in the same basis by 
projecting the self energy in the localized Wannier ba- 
sis |xm(k)), S mm , to the basis of Bloch states \ip u (k)), 
IW'(k), 



(1) 



where P mi/ (k) = (x m (k)|Vv(k)). From the diagonal ele- 
ments Yi ul , we obtain the mass enhancements of the re- 
spective Fermi surface pocket at the fc-points indicated 
in Fig. [3j the values are given in Table [TT] 



TABLE II. Mass enhancements of the maximal/minimal 
Fermi surface orbits in LiFeP. The mass enhancements are 
measured at the fc-points indicated in Fig. [3] 



orbit 


la/b 


2a/b 


3a/b 


4a/b 


5a/b 


pocket 


inner 


middle 


outer 


outer 


inner 




hole 


hole 


hole 


electron 


electron 


m* 
m LDA 


1.47/1.34 


1.48/1.47 


1.69/1.70 


1.49/1.46 


1.52/1.37 



The significantly smaller mass enhancements of the 
middle hole pocket (orbit 2) measured in the quantum os- 
cillation experiments are not seen in LDA+DMFT. This 
suggests that this pocket is differently affected by the 
coupling to some scattering channel like spin fluctuations 
or nonlocal correlations which are not captured by our 
LDA+DMFT approach. 

In summary, we reported LDA+DMFT calculations on 
LiFeP and LaFePO where we find a change of the Fermi 
surface topology upon inclusion of correlations in both 
compounds, namely the opening of an outer hole pocket 
at r in LaFePO and the opening of an inner hole pocket 
at Z in LiFeP, both with Fe d z 2 orbital character. As dis- 
cussed by Kemper et alP, this might promote the nodal 
gap and weak pairing strength, i.e. low T c , in these mate- 
rials. Whereas this pocket has been observed in ARPES 8 
for LaFePO, the current experimental situation for LiFeP 
doesn't allow for a conclusive testing of our predictions 
and further experimental work is desired. Also we find 
that the pecularities of the middle hole pocket in LiFeP 
observed in dHvA experiments and not reproduced in our 
LDA+DMFT approach reveal the importance of scatter- 
ing channels beyond local correlations. 

We would like to thank A. Coldea, C. Putzke, A. Car- 
rington, and P. J. Hirschfeld for useful discussions and 
the Deutsche Forschungsgemeinschaft for financial sup- 
port through grant SPP 1458, the Helmholtz Association 
for support through HA216/EMMI and the centre for sci- 
entific computing (CSC,LOEWE-CSC) in Frankfurt for 
computing facilities. 



5 



* ferber@itp.uni-frankfurt.de 

1 A. F. Kemper, T. A. Maier, S. Graser, H.-P. Cheng, P. J. 
Hirschfeld, and D. J. Scalapino, New J. Phys. 12, 073030 
(2010). 

2 G. Wu, Y. L. Xie, H. Chen, M. Zhong, R. H. Liu, B. C. 
Shi, Q. J. Li, X. F. Wang, T. Wu, Y. J. Yan, J. J. Ying, 
and X. H. Chen, Journal of Physics: Condensed Matter 
21, 142203 (2009). 

3 Y. Kamihara, H. Hiramatsu, M. Hirano, R. Kawamura, 
H. Yanagi, T. Kamiya, and H. Hosono, J. Am. Chem. 
Soc. 128, 10012 (2006). 

4 T. M. McQueen, M. Regulacio, A. J. Williams, Q. Huang, 
J. W. Lynn, Y. S. Hor, D. V. West, M. A. Green, and 
R. J. Cava, Phys. Rev. B 78, 024521 (2008). 

5 M. Yamashita, N. Nakata, Y. Senshu, S. Tonegawa, 
K. Ikada, K. Hashimoto, H. Sugawara, T. Shibauchi, and 
Y. Matsuda, Phys. Rev. B 80, 220509 (2009). 

6 J. D. Fletcher, A. Serafin, L. Malone, J. G. Analytis, J.- 
H. Chu, A. S. Erickson, I. R. Fisher, and A. Carrington, 
Phys. Rev. Lett. 102, 147001 (2009). 

7 C. W. Hicks, T. M. Lippman, M. E. Huber, J. G. Analytis, 
J.-H. Chu, A. S. Erickson, I. R. Fisher, and K. A. Moler, 
Phys. Rev. Lett. 103, 127003 (2009). 

8 D. H. Lu, M. Yi, S.-K. Mo, A. S. Erickson, J. Analytis, 
J.-H. Chu, D. J. Singh, Z. Hussain, T. H. Geballe, I. R. 
Fisher, and Z.-X. Shen, Nature 455, 81 (2008). 

9 M. M. Qazilbash, J. J. Hamlin, R. E. Baumbach, L. Zhang, 
D. J. Singh, M. B. Maple, and D. N. Basov, Nat. Phys. 5, 
647 (2009). 

A. I. Coldea, J. D. Fletcher, A. Carrington, J. G. Analytis, 
A. F. Bangura, J.-H. Chu, A. S. Erickson, I. R. Fisher, 
N. E. Hussey, and R. D. McDonald, Phys. Rev. Lett. 101, 
216402 (2008). 

S. Suzuki, S. Miyasaka, S. Tajima, T. Kida, and M. Hagi- 
wara, J. Phys. Soc. Jpn. 78, 114712 (2009). 

2 S. L. Skornyakov, N. A. Skorikov, A. V. Lukoyanov, A. O. 
Shorikov, and V. I. Anisimov, Phys. Rev. B 81, 174522 
(2010). 

3 Z. P. Yin, K. Haule, and G. Kotliar, Nature Mater. 10, 
932 (2011). 

4 J. H. Tapp, Z. Tang, B. Lv, K. Sasmal, B. Lorenz, P. C. W. 
Chu, and A. M. Guloy, Phys. Rev. B 78, 060505 (2008). 

5 Z. Deng, X. C. Wang, Q. Q. Liu, S. J. Zhang, Y. X. Lv, 
J. L. Zhu, R. C. Yu, and C. Q. Jin, Europhys. Lett. 87, 
37004 (2009). 

6 K. Hashimoto, S. Kasahara, R. Katsumata, Y. Mizukami, 



M. Yamashita, H. Ikeda, T. Terashima, A. Carrington, 
Y. Matsuda, and T. Shibauchi, Phys. Rev. Lett. 108, 
047003 (2012). 

17 K. Kuroki, H. Usui, S. Onari, R. Arita, and H. Aoki, Phys. 
Rev. B 79, 224511 (2009). 

18 H. Ikeda, R. Arita, and J. Kunes, Phys. Rev. B 81, 054502 

(2010) . 

19 R. Thomale, C. Piatt, W. Hanke, and B. A. Bernevig, 
Phys. Rev. Lett. 106, 187003 (2011). 

20 For a review see P. J. Hirschfeld, M. M. Korshunov, and 
I. I. Mazin, Rep. Prog. Phys. 74, 124508 (2011). 

21 S. Kasahara, K. Hashimoto, H. Ikeda, T. Terashima, 
Y. Matsuda, and T. Shibauchi, Phys. Rev. B 85, 060503 
(2012). 

22 C. Putzke, A. I. Coldea, I. Guillamon, D. Vignolles, A. Mc- 
Collam, D. LeBoeuf, M. D. Watson, I. I. Mazin, S. Kasa- 
hara, T. Terashima, T. Shibauchi, Y. Matsuda, and 
A. Carrington, Phys. Rev. Lett. 108, 047002 (2012). 

23 M. Aichhorn, L. Pourovskii, and A. Georges, Phys. Rev. 
B 84, 054529 (2011). 

24 P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and 
J. Luitz, WIENZk, An Augmented Plane Wave + Local Or- 
bitals Program for Calculating Crystal Properties (Techn. 
Universitat Wien, Austria, 2001). 

25 P. Werner, A. Comanac, L. de' Medici, M. Troyer, and 
A. J. Millis, Phys. Rev. Lett. 97, 076405 (2006). 

26 B. Bauer, L. D. Carr, H. G. Evertz, A. Feiguin, J. Freire, 
S. Fuchs, L. Gamper, J. Gukelberger, E. Gull, S. Guertler, 
A. Hehn, R. Igarashi, S. V. Isakov, D. Koop, P. N. Ma, 
P. Mates, H. Matsuo, O. Parcollet, G. Pawlowski, J. D. 
Picon, L. Pollet, E. Santos, V. W. Scarola, U. Schollwock, 
C. Silva, B. Surer, S. Todo, S. Trebst, M. Troyer, M. L. 
Wall, P. Werner, and S. Wessel, J. Stat. Mech. , P05001 

(2011) . 

27 E. Gull, P. Werner, S. Fuchs, B. Surer, T. Pruschke, and 
M. Troyer, Comp. Phys. Commun. 182, 1078 (2011). 

28 V. I. Anisimov, F. Aryasetiawan, and A. I. Lichtenstein, 
J. Phys.: Condens. Matter 9, 767 (1997). 

29 M. T. Czyzyk and G. A. Sawatzky, Phys. Rev. B 49, 14211 
(1994). 

30 L. de' Medici, Phys. Rev. B 83, 205112 (2011). 

31 S. V. Borisenko, V. B. Zabolotnyy, D. V. Evtushinsky, 
T. K. Kim, I. V. Morozov, A. N. Yaresko, A. A. Kordyuk, 
G. Behr, A. Vasiliev, R. Follath, and B. Biichner, Phys. 
Rev. Lett. 105, 067002 (2010). 

32 C. Putzke, private communication. 



